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in the experiment ten sets of 50 scattered pulses of pulse duration 88 pis and 180 
kHz carrier frequency were averaged, and the envelope of the decaying tail was fitted 
with a logarithmic curve. The decay constants for the various sediments showed a 
variation with the grain size expressed by a = 10.3 - n/0.0035 where a is the grain 
diameter in mm and yn the decay constant in Np/pts. The value of a can be estimated 
within + 2.8 mm for a decay constant of 2.5 x 1072 Np/pts. This relation is limited to 
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I. INTRODUCTION 


The first use of sound to detect underwater objects by echo-ranging was in World 
War I, but it was not until World War II that sonar development became significant. 
Detection and ranging are the main functions of sonar systems, which include the 
fathometer, a widely used and very important application. The fathometer measures 
the round-trip travel time of short sound pulses between a transmitter located at the 
Waters surface and the bottom. If the sound speed in the water is known, depth 
measurements are possible and the safety of navigation in near-shore areas is 
improved. 

The fundamental principle is simple but the process involves some complexity 
due to the variation in sound speed and attenuation in the water. These processes are 
inherent to the water, but when a target or the bottom 1s present, reflection, absorption 
and incoherent scattering at the interface between the target or bottom must be 
considered. 

This is the situation present in a hydrographic survev operation: A transducer in 
the water sends sound pulses to the bottom which are reflected back. Along the sound 
path transmission losses due to spreading, absorption and scattering occur. Changes in 
temperature, salinity and pressure along the path of a sound rav cause refraction, 
Which is not a very important problem in hydrography, because the transmitter points 
downward. Upon reaching the bottom, sound is both reflected and transmitted into 
the bottom so that its amplitude is decreased. Some of the sound transmitted into the 
bottom is scattered backward where it combines with the reflected pulse, distorting the 
initial pulse shape. This pulse distortion is the subject of this paper. 

In hydrography it is the front edge of the received pulse which is normally of 
interest, since it gives the distance to the bottom. In addition, however, more 
information may be obtained from the same received pulse. In fact, the tail of the 
received signal can give information about the nature of the bottom. There is thus the 
possibility to remotely classify the bottom type, an operation essential to the 
preparation of nautical charts. 

The intent of this project is to verify a laboratory technique previously studied 
and established by Diaz (1986) and Chang (1986), in which the trailing edge of the 


1] 


reflected pulse is characterized in terms of a Straight line fitted to the envelope by a 
least square method. In this study the same experimental procedures are applied to 
natural sediments and artificial materials to test the theory relating the tvpe of bottom 


sediment to the received echo level. 


If. THEORETICAL CONSIDERATIONS 


A. DEFINITION OF REVERBERATION 

In any echo-ranging or listening device the desired signal is always received in the 
midst of a certain amount of extraneous noise (National Defense Research 
Commutee,1969). Particularly in underwater acoustics, this noise can arise from 
different sources such as breaking waves, snapping shrimp, surf and shipping, which 
combine to produce broadband ambient noise (Kinsler, Frey, Coppens and Sanders, 
1983). 

However if dealing with an active system like the sounding equipment used in 
hydrography, another type of interference arises, which is directly due to the 
transmission of the pulses into the water. This component is called reverberation and 
results from the existence of particles or discontinuities in the path of the pulse that act 
as small scattering centers, reradiating the acoustic energy and returning a weak echo 
to the transducer. The combination of all these small echoes makes the reverberation. 
Like ambient noise, reverberation is usually undesired, because it masks the reflected 
signal. 

1. Surface, Bottom and Volume Reverberation 
When small scatterers or discontinuities are observed along a reflecting 
surface, like the water-air interface, surface reverberation occurs. It is important when 
a reflecting surface has a certain roughness compared to the wavelength of the sound. 
When the reflection is produced at the bottom surface, and if this bottom surface 
shows some irregularities or small objects, then bottom reverberation is produced. But, 
if the reverberation is due to scatterers existing in the body of the propagation medium, 
i.e. the water column, then we have volume reverberation. In this project the volume 
reverberation produced in the sediment is of interest because, as it was pointed out 
before, the amount of volume reverberation observed can give an indication of the type 
of sediment present. This suggests the possibility of rapidly determining the types of 
bottom sediments present during hydrographic soundings. 
2. Reverberation Level 
If an active system with an intensity Ne at the transmitter is considered, with a 


target at range r, the intensity level, IL(r), of the signal reaching the target 1s 


[L(r) = 10 log [I(r)/I_.4q dB (250) 
where 
oS (2.2) 


for transmission loss by spherical spreading. 

After being scattered by a volume, dV, the intensity level depends on the ability of the 
target to reflect the incident energy. This ability is expressed by S, the scattering 
strength. | 


ST Oras (2.3) 
where s, is the scattering coefficient per unit volume 
Se 0/47 (2.4) 


and o 1s the backscattering cross section of the unit volume. Thus, the scattered 


intensity, I(r) at 1 m from the target is 
I(1) = (Ipir2)s,V (2.5) 
S O° Pv 


where V is the volume occupied by scatterers that contribute to the reverberation 
received at the receiver at time t. One essential step in obtaining the reverberation 
level (RL) is to compute the scattering volume which is a function of the transmitting 
pulse length (t), and the directivity of the transducer (here assumed monostatic). The 


geometry occurring in hydrographic surveying is shown in Figure 2.1. 





Figure 2.1 Echo sounding geometry and elementary portion 
of volume reverberation contribution. 


The thickness of the reverberating volume must be such that all the scatterers in it 
contribute reverberation at the same time at the receiver, (Kinsler, Frey, Coppens and 


Sanders, 1983) and that depends directly on the pulse length or 
L. =-cu2 (2.6) 


where Tt is the pulse duration. For a geometric cross section of the beam pattern, A, 


we have 
V = Ac(t/2) (2.7) 
Substituting this into Equation 2.5 gives 
1.(1) = (Ip/r)s ,Ac(t/2) (2.8) 
and the intensity at the receiver becomes 
I. = (Ip/r4)s,Ac(t/2) (2.9) , 


The volume reverberation level in the form of a sonar equation is 


RL = SL- 40 logr + S, + 10 log V (2.10) 


Urick (1983) suggests that in reverberation backgrounds it is easy to find the echo-to- 


reverberation ratio for a target of known target strength: 


EL = SL-2TL + TS (2.11) 
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where 

EL = echo level 

SL = source level 
TS = target strength 


Subtracting Equation 2.10 from Equation 2.11 we obtain 


ei = TS = (S. ae ocay) (23) 


but (S,. + 10 log V) is the reverberation target strength, and the echo-to-reverberation 
ratio is only the difference of two strengths, one for the reflected echo and the other for 
Ememscattcred One. 

How can this suggestion be applied to our project ? In it we are interested in 
finding the RLs produced by different materials. Considering the reflection only at the 
fluid-sediment interface (no incoherent surface or volume scattering), the reflected 
signal will have the same shape, no matter which material is below that interface, being 
reduced in amplitude onlv by the reflection coefficient of the sediment. 

If the incoherent scattering produced by the sediment particles is added the 
received scattered signals will differ in shape, (mainly in the trailing edges) for different 
sediments. These differences are due to the incoherent scattering represented in the 
sonar equation by RL. 

There is no way to measure directly the pure (coherent) pulse reflected at the 
water-sediment interface because the reflected signal is “contaminated” by incoherent 
scattering, but it can be approximated by the pulse reflected from the water-air 
interface. The only difference between these is the voltage amplitude because the 
specific acoustic impedances are different, but this can be corrected by normalizing the 
received voltages. 

In Chang (1986) this pure, reflected signal was obtained using a transducer 
lying on the bottom and transmitting upward to the water-air interface. The signal 
received is then purely reflected, assuming that no surface reverberation is produced at 
the interface or that its magnitude is negligible compared to the volume reverberation 
produced in the sediment. This assumption also applies to the water-sediment 


interface, which must be perfectly level and smooth. 
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By using the same transmitting conditions, but with the transducer located at 
the water-air interface and transmitting downward to the bottom, a different signal is 
received back in the transducer; this represents the volume reverberation. 

Geometrically this procedure can be illustrated with the pictures of the two 
distinct waveform envelopes. Besides the difference in amplitude, the most noticeable 
feature is that the received waveform shows a slower decaying of its tail, where 
reverberation is geometrically expressed. Mathematically, the decaying tail can be 


represented by 


V(t) = Vo e*' (2.13) 


for the water-air interface and 


V(t) = VoeM (2.14) 


for the water-sediment interface, where ve is normalized to the same value in each case 
(Chang,1986). The difference between these two expressions gives a measure of the 


amount of volume reverberation coming from that particular sediment. 


AV(t) = V{(e NM - e%) (255) 
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I‘igure 2.2 Geometric interpretation of volume reverberation. 


If. BACKGROUND 


Many research projects have been devoted to the study of reverberation in the 
ocean. However, most have been concerned with scattering produced in the water 


column and only a few with volume reverberation in marine sediments. 


A. LITERATURE REVIEW 

Urick (1953) studied backscattering from the bottom as a function of pulse 
length, frequency and grazing angle at several locations in a harbor using a tiltable 
transducer mounted on a barge. The results were expressed in terms of scattering 
strength per unit area of bottom. Mackenzie (1961) made bottom-reverberation 
measurements using two low frequencies and different pulse lengths as well as different 
grazing angles. He analysed the data to obtain scattering coefficient for the bottom. 
McKinney and Anderson (1963) made the same measurements but for three different 
types of sediment: mud, sand and gravel. They concluded that the bottom sediment is 
a major factor in the backscattering of sound and that the knowledge of particle size 
distribution of the sediment is useful in estimating the reverberation level. No model 
was formulated. 

Recent relevant work (Dodds [{1982,1984], Breslau [1967] and Meng [1982]) deals 
with the problem of using backscattering measurements to remotely classify bottom 
types. Dodds used a broadband acoustic source (1-10 Khz) directed vertically to the 
sea floor and obtains sonograms from which volume scattering and surface roughness 
parameters are obtained. Meng (1982) discussed the classification of seabed sediments, 
based on the echo envelope of normally reflected high-frequency sound pulses. This 1s 
apparently similar to the present work. It would be interesting to investigate his 
procedure, processing and achievements. However this paper was published in 
Chinese, and only a short abstract is available in English. 


B. VOLUME REVERBERATION MODELS 

The backscattering of sound from the ocean bottom has been studied for many 
years, but apparently only Morse and Ingard (1968) and Urick (1983) have developed 
theoretical models appliable to the present study. 
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These models have, however, a restricted application since they make several 
assumptions that are not satisfied in the present case, where the wavelength is 8 mm 
and the average diameter of the scatterers ts between 0.3 and 5.3 mm. 

Morse and Ingard’s assumptions and their relevance to the present study are: 

1) Scattering objects are small compared with the wavelength. 
(satisfied here onlv for the fine sand). 

2) Dimensions of the scattering region are much larger than the 
wavelength. 

(Assumption satisfied for all four sediments). 

3) The scatterers are all spheres. (Not too bad an assumption for all 
four sediments). 

4) The scatterers are sparsely packed and no multiple scattering is 


considered. (Very poor assumption for all sediments). 


For monostatic sonars Morse and Ingard give the volume scattering coefficient 


per unit volume as 
ne (2m) Ko ay ia (3.1) 


where “a” is the radius of the scatterers, N the number of scatterers per unit volume 


and k is the wave number, 

Ne = (Ky-Ko) Ko 

Where K,, is the compressibility of the scatterer, K,, the compressibility of the fluid and 
1p = GPn-3Po) / (2Pn-Po) 

where p,, is the scatterer density and p,, the fluid density. 


The scattering coefficient per unit volume as given by Urick (1983) 1s 
s,=No (2-2) 
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where 6, the backscattering cross section, is given by 
o = 2.8(ma2\(k a) (3.3) 


where “a” is the scatterer dimension. This model is good for a single sphere (no 


multiple scattering) and for values of Ka < < 1. 


C. DERIVED MODEL 
The present project is an extension of the recent work developed by Diaz (1986) 
and Chang (1986). Chang measured volume reverberation on aggregate and fine sand 
to compute the corresponding volume-scattering coefficients. Applying Morse and 
Ingards (1968) and Urick’s (1983) models Chang did not get satisfactory results, 
verifying that these models are not appropriate for ocean sediments. A simple model 
was then set up: 
The envelope of the electrical pulse shape can be described by a square wave of 


pulse duration T 


I(t) = TG a(t) (3.4) 


where 
I(t) = intensity at time t after the beginning of the pulse 
I, = constant intensity amplitude 

t = pulse length 


H(t) = unit step function shown in Figure 3.1 a) 


Ost) 
H(t) = 
1 t>0 


H(t-t) is then: (Figure 3.1 b)) 


0 t-T<O 
H(t-t) = 
1 t-t>Q 


22 


The transmitted pulse is not exactly square due to damping of the transducer, but it 
is tranformed into an acoustic pulse with an exponential rise, a flattened top and an 


exponential decay. This shape can be described by: Figure 3.1 c¢ 


Rt) = (ee @*y (ry + (1-e°@9) H(t-r) (3.5) 





Figure 3.1 Step unit function and transmitted acoustic pulse. 


Where t is the time measured from the beginning of the pulse and a@ the rise and decay 
constants. For a transducer located in the water sending pulses towards a smoothed 
surface the intensity received by the same transducer must be according to Chang 
(1986) 


P= (fT 2 sy ePS 1, e-2b ley (r+ 8)? + Rly [repr (3.6) 


Wwiere: 


T; = intensity transmission coefficient 
B = absortion cofficient in the sediment = 0.25 f (Urick,1979) 
& = affected sediment depth = c t/2 
R; = intensity reflection coefficient 
r = distance from transducer to the sediment surface 
The first term inside the square brackets represents the scattered component of the 


scattered intensity for the volume of the sand and the second the reflected one. 


Forr>>& 

Le = (Ty? sy Tyr?) ferPS r-28/c)ae (3.7) 
where 

flt-2& ic) = [Leet 25/C yp (1.28/c) + [Lae AEE 25/ Cy (1-2-26/c) (3.8) 


After integration between & = 0 and & = ct/2 


[, = A [(1-2°CBU2) p.(e-CBUi2. et) B + 2a/c) (3.9) 


te (e7CB(t-T)/2_ 1/8 ae (e7CB(t-t)/2_,-a(t-T)) / (-B + 2a/c)] 
where 


are 4 
A=TisyIo/t (3.10) 


With some algebraic operations 
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I, = A(1/Pfe"Bt/2(ecBt/2. 1 (3.11) 


ae 1/(-B orks 2a/c)fe"CPt!2(ecBt/2_4) + enun@ieens “)} 


and the reflected component 


a die: 
P= 1,Rit (ei9*-1) e@! (3.12) 


Adding the two components incoherently, the total (reflected and 


scattered) intensity received by the transducer will be 


[= B{T,2s, /Ple CPU 2(eCBt/2. 1) 4 T.4s,/(-B+ 2a/c) (3.13) 


fe"CBt/2/QcBt/2. 4) 4 en Fly AT) 4 Ri(e° AT 1)e° AY 


where B = ae 
The expression looks quite complex but all the variables are known except s,. If 


the part included inside braces is called C, then 
cc 


The quantity inside brackets is the received intensity considering only spherical 
Spreading and C accounts for the scattered and reflected intensity. The only interesting 
part is for t > 7, that corresponds to the decay portion. That is also what is measured 
in the experiment, i.e. the decay of the pulse tail. 

To find the value of Sy Chang combined the expression for C with his 


the corresponding values for C were 
2 
C ' 


experimental results. Guessing values for s, 
obtained. He repeated then the arithmetic procedure until the value of the log 
that expresses the slope decay became close to the decay constant obtained 
experimentally. The value guessed for s,, is then the volume scattering coefficient for 


the sediment studied. 
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This simple model appears to approximate the volume scattering. Further study 


of volume reverberation from ocean sediments to improve the mode 


l was thus 
suggested. 
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IV. EXPERIMENT 


A. INSTRUMENTATION 
1. Experiment Environment 

All of the preliminary and main experiments were conducted in two steel- 
bound glass tanks filled with fresh water and having interior measurements of 0.7 m by 
0.7 m by 0.6 m. Four different materials, fine sand, aquarium sand, aggregate 

and glass beads, were placed at the bottom of the tank to serve as “sediments”. 

The tank had 0.35 m of water overlying a 0.25 m thick layer of “sediment”. 
Specifically, the aggregate and the glass beads were kept inside a plastic bucket during 
the experiment to reduce the amount of material used. Because the transducer beam is 
narrow (projected spot is a circle with 0.1 m radius at a 0.3-m range), the beam fell 
completely inside the bucket (Figure 4.1). 

The fine sand, also studied by Bradshaw (1981), Diaz (1986) and Chang 
(1986), had an average size of 0.3 mm and varied from 0.15 mm to 0.7 mm. The 
aggregate, also used by Diaz (1986) and Chang (1986), was less homogeneous in size 
and composition than the fine sand and had an average size of 5.3 mm. 

It was convenient to make measurements using a sediment size between 0.3 
mm and 5.3 mm for a better evaluation of the volume reverberation behavior with the 
grain sediment size, and for that a third sediment was selected, the aquarium sand. This 
one was sieved and sized to an average size of 1.72 mm and varied from 0.84 mm to 
3.36 mm with many inhomogeneities in composition and shape. It was expected that 
these variations in size, shape and constitution for the same sediment would raise some 
difficulties and this was confirmed during the experiment. The idea of performing the 
experiment on an ideal homogeneous sediment lead to the use of spherical glass beads, 
whose density (2.5 g/cm?) is not far from the real sediments’ densities and size is 1.7 + 
0.3 mm. 

The four different materials, namely the fine sand, the aggregate, the aquarium 
sands and the glass beads, were kept in water for two to three days before 
measurements were made. A water jet from a hydraulic pump, or manually stirring 
was used to remove trapped air bubbles before and during the measurement process. 
Bleach was added to the water to control biological growth, and a wetting agent was 


also mixed with water to help remove air bubbles. 
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Figure 4.1] 


Experimental Geometry Design. 


2. Electronic Equipment 

A General Radio Model 1310 Oscillator generated a 181-kHz signal that was 
simultaneously fed into a frequency counter and a General Radio type 1396-A Tone 
Burst Generator to produce a 16-cycle pulse, which was then amplified by a Hewlett 
Packard 467A power amplifier before passing to a Datasonics Transmit) Receive (T/R) 
switch that controlled the transmission and reception by the same transducer. The 
transmitted signal was used to trigger an analog oscilloscope (Tektronix Type RM 503) 
and a Nicolet model 3091 Digital oscilloscope. 

The transducer, a Naval Research Laboratory F-41 circular piston type with 
an $.8-cm diameter active face, has a resonance frequency of 182 kHz and a half-beam 
width of 10 degrees (Diaz,1986). 

The signal received by the transducer was amplified 20 dB by a 465A Hewlett- 
Packard pre-amplifier and then filtered by a Spencer Kennedy Laboratories, Inc. Model 
302 electronic filter, used as a high pass filter to avoid low frequency noises. The 
filtered signal was finaliy passed to the oscilloscope which was controlled by a 3421A 
Hewlett Packard Data Acquisition/Control unit attached to an Hewlett Packard 86 
desk-top computer. 

All the equipment worked reasonably well : The frequency generator showed a 
drift of + 150 Hz during a complete set of 10 measurements, but this did cause any 
significant changes in the pulse waveform, as long as the transducer was operated near 
its resonance. The digital oscilloscope was controlled by the computer and the the 
data acquisition/control unit, but sometimes during the measuring process the 
oscilloscope stayed locked in the STORE mode making it necessary to interrupt the 
measurements and reset it. 

Each set of measurements was the result of 50 samples and took about 45 
minutes. A faster sampling rate could speed up the measuring process and probably 
give us a more accurate waveform reconstitution. A block diagram of the equipment 1s 


shown in Figure 4.1. 


B. PRELIMINARY MEASUREMENTS 
Preliminary measurements of density, speed of sound and reflection coefficient 


were made for later computations of acoustic characteristics. 
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Figure 4.2 Electronic Equipment block diagram. 





1. Density 
Bradshaw (1981) measured the water-saturated density, 
sand and obtained: 
Paix = 1.98 + 0.03 g/cm! 
p = 2.69 + 0.01 gicm? 


For the aggregate Diaz (1986) obtained : 
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Pix» and dry density, p, of fine 


P mix = 1.97 g/cm 
p = 2.64 g'cm? 
Density measurements were also performed for the aguariumsand and for the glass 


beads. The results are summarized in Table }. 


2. Speed of Sound 

The measurement of sound speed in the sediment was performed inside the 
glass tank using three LC10 transducers, one acting as a transmitter and the other two 
as receivers. The transducers were buried in the sediment at the same depth. The 
transmitter and the closest receiver were kept fixed and the third transducer was moved 
to different positions along the line connecting the first two transducers. 

The difference in reception time at the two receivers was measured by means 
of the digital oscilloscope (Figure 4.3), and by measuring the distance between the 
receivers it was possible to calculate the sound speed. Fifteen measurements were 
made on the aquarium sand resulting ince = 1588 + 22 m/s. 

Hamilton et al. (1956) studied the physical properties of marine sediments 
from off the coast of San Diego, California, and established some relations between the 
sediment grain sizes and their physical properties such as porosity, density, speed of 
sound, etc. They observed a general increase of sound speed with increase in medium 
grain size. Thus, the value of the sound speed for the aquarium sand was smaller than 
the 1610 m/s obtained by Bradshaw (1981) for the fine sand, and does not agree with 
Hamilton et al.. the last conclusion. To investigate the consistency of that value, a 
second set of measurements was made and the result was c = 1590 = 42 ms, 
confirming the first one. 

The same experimental procedure was conducted on the glass beads and the 


results were c = 1822 = 48 m/s 
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Figure 4.3 Detail of speed of sound measurement procedure. 


3. Reflection Coefficient 
The reflection coefhicient for sound rays passing from one medium to another 
is the ratio between the reflected and the incident acoustic pressures, and for normal 
incidence it depends on the specific acoustic impedance pc, of each medium. The 
reflection coefficients for aquarium sand and glass beuds were measured as the ratio 
between the maximum voltage amphtude in the signal reflected from the sediment and 


the signal reflected by the water-air interface at the same distance. 


R = VV; (4.1) 


The reflection coefficient can also be computed using 


R = (r5 = ry dry a5 ry) (4.2) 





where 


ry} = Pc, is the specific acoustic impedance for the water 


I> = Pc for the sediment 


The measured, ae and computed, Sex values of the reflection coefficient, as 


well as other physical characteristics are summarized in Table 1. 


TABLE 1 
SEDIMENT PHYSICAL CHARACTERISTICS 





fine 
sand 


aggregate 


aquarium 
san 


lass 
wena 


| a 0.3 5.3 1.72 Mi 
(mm) (+ 0.3) (+ 0.51) (+ 0.24) (+ 0.3) 
9 2.69 2.64 2.59 2.50 
(g’cm?) (+ 0.01) (+ 0.05) (+ 0.05) (+ 0.05) 
a 1.98 1.97 2.02 1.91 
(¢, cm?) (+ 0.03) (+ 0.05) (+ 0.02) (+ 0.01) | 
¢ 1610 1555 1588 1822 
(my So) (= 22) (+ 48) 
Rie 0.36 0.36 0.36 0.53 : 
(+ 0.09) (+ 0.05) | 
R. 0.38 0.52 0.38 0.41 | 
(+ 0.01) (+ 0.01) 
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C. MEASUREMENT PROCEDURES 
1. Experiment Technique 

To measure the exponential decay of the tails of the pulses scattered from the 
sediment the experimental procedure of Diaz (1986) and Chang (1986) was used. The 
transducer was clamped to a horizontal wooden rail that kept the transducer level with 
its active face pointing downward to the sediment as shown in Figure 4.1. 

In his experiment Chang observed large variations in the received pulse shapes 
between different samples of the same sediment and concluded that the variations were 
due to inhomogeneities in the sediment. Because of this, he did not obtain consistent 
results. To smooth the variations in the results a set of ten measurements was 
performed with various transducer positions, and the results were averaged 

For the aggregate sand it is not difficult to accept that important 
inhomogeneities are present in its components, but the /fime sand, looks quite 
homogeneous. However, the waveforms did show large variations. By performing the 
experiment using a thin metal plate covering the sediment, Chang showed that those 
inhomogeneities were interior to the sediment. 

To investigate these facts further, the first experiment performed in this 
present work consisted of keeping the transducer in a fixed position and, between each 
measurement the sediment was stirred using a water jet produced by a pump. During 
this stirring process clouds of air bubbles were released from within the sediment, and 
after that, the sediment texture changed completely, becoming much “silkier” to the 
touch. 

The reflected signals then became consistent, as can be seen in Figure 4.4, where pulses 
scattered from this sediment, after each stirring, are displayed. The inhomogeneities 
found by Chang in his experiment appear to have been caused by occluded gas within 
the sediment. Thereafter, each sample of sediment being studied was. stirred 
periodically, to remove any air bubbles. 

2. Measurement Considerations 

The measurements made on the aquarium sand with this new tecnique still 
revealed some variation, mostly in the last part of the pulse tails, as can be seen in the 
group of five waveforms in Figure 4.5 and in the measured decay slopes presented in 
Ghaptens: 

Measurements for the aggregate were repeated to determine if the 


inconsistencies in our results were also due to the presence of gas bubbles within the 
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Figure 4.5 Waveforms of pulses scattered from aquarium sand. 
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Figure 4.6 Waveforms of pulses scattered from aggrcgate. 
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saniple or due to the visible physical inhomogeneities present. This time, the variations 
remained (Figure 4.6) although when the sediment was deeply stirred only a few 
bubbles were observed during the stirring. A conclusion can be drawn from these 
experiments with the different sediments : The variations in the tail shape of the pulse 
increases as the grain size of the sediment increases. 

The consistent dependence of the results on the grain size, indicated that we 
may have been working with grain dimensions too large relative to the wavelength. In 
fact the average dimension of the aggregate is 5.3 mm and the wavelength is 8 mm, 
values that are very close. One of the assumptions necessary to determine the volume 
reverberation is that the surface scattering is negligible, because the sediment surface 1s 
levelled and smoothed. The levelling of the bottom surface was made by using a device 
that, rolling on the tops of the tank walls, produced a horizontally levelled bottom 
surface. Even with levelling, when the grain sizes are relatively big there is always a 
gap between adjacent grains that destroy the smoothness and produces the effect of a 
rough surface. 

A simple experiment was performed to determine the importance of surface 
scattering. A Styrofoam plate, to simulate a pressure release interface was weiglhited 
down on the sediment surface and covered by a thin laver of sediment, and the 
scattered pulses were analysed. In this case, anv changes in the pulse waveforms would 


result, not from volume reverberation, but from surface scattering. 
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Figure 4.7 Waveform of pulse reflected from a thin laver of aggregate. 


It was sufficient to observe the waveform on the oscilloscope to find that the 
decay changed drastically each time the transducer was slightly moved. Merely moving 
some grains from the thin layer caused the received pulse to change completely. The 
pulses (Figure 4.7 to Figure 4.9) showed many of the same features as were observed 
during the volume reverberation measurements. It was concluded that surface 
scattering was a significant factor in reflection from sediment when the grain size was 
comparable to the wavelength of the sound. 

Performing the same experiment with Aquarium sand led to a similar 
conclusion that surface scattering was playing an important role. Although the pulse 
waveforms did not show large variations, the values for the decay slopes showed a 
significant spread, indicating the presence of surface scattering. 

To analyse the extent of the influence of inhomogeneities on the results 


obtained in the measurements performed with the three previously indicated sediments, 
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Figure 4.8 Waveform of pulse reflected from a thin layer of aggregate. 


another set of ten observational runs was made using glass beads. Jt was expected that 
the reflected pulses could be more consistent due to their homogeneity. This did not 
happen, however, and Figure 4.11 shows five different waveforms reflected by the glass 
beads. They present approximately the same features as the pulses reflected by the 
aquarium sand whose average grain size is close to that of the beads dimensions, 
showing once more the influence of the grain size in its scattering properties. 
3. Data Processing 

The computer programs used in this project were designed and tested by 
Chang (1986). They were used to obtain the average slope for each sample: 
THESTS3 ts a program designed to convert and normalize the waveform data obtained 
from the Nicholet digital oscilloscope. Since the pulses contain positive and negative 
values, these are squared and a series of fifty waveform samples are averaged to give an 


accurate representation of the waveform envelope. Because of jitter in the trigger 
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Figure 4.9 Waveform of pulse reflected from a thin layer of aggregate. 


signal, each “bin” of the scope will sample the waveform at various times during a 
cvcle, and the resulting curve will be a representation of the envelope of the pulse. 

The result can be seen in Figures 4.12 to 4.14 where four tails of the puises 
reflected from the three different sediments are plotted. 
MAXIMUM is a program used to choose the local maxima for the curves obtained 
from THESIS3 stored on disk. The square root of these maxima is then taken and 
stored in the same disk. The result is shown in Figures 4.15 to 4.17 which show 
respectively the envelopes of the pulse tails shown in Figure 4.12 to Figure 4.14. 
PLOT3 is used to produce the plots of the waveform envelopes and to apply a least 
square method to find the best fitted straight line, its slope and correlation coefficient. 
AVERAGE. The measurement of the slope of the tail was very sensitive to any 
change in the sediment sample, and to get a consistent result ten sets of ineasurements 
were perforrned, and this program is applied to compute a chain average in the 


following wav: 
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Figure 4.10 Waveform of pulse reflected from Styrofoam plate. 
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Figure 4.11 Waveform of pulses scattered by glass beads. 
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Figure 4.12 Tail waveform of pulse reflected from fine sand. 
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Figure 4.13 Tail waveform of pulse reflected from aquarium. 
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Figure 4.14 Tail waveform of pulse reflected from aggregate. 
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Figure 4.15 ‘ail envelope of a pulse reflected from fine sand. 
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Figure 4.16 Tail envelope of a pulse reflected from aquarium. 
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Figure 4.17 Tail envelope of a pulse reflected from aggregate. 
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Each averaging step is divided by the number of terms envolved so that each 
data set has the same weight. The DAT10 is then the file that contains the averaged 
maxima that plotted in Figure 4.18 shows a much more regular envelope than the 


envelopes of the local maxima plotted, for example, in Figure 4.12. 
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Figure 4.18 Final envelope of an averaged pulse. 
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V. RESULTS 


The decay constants of the trailing edge of sound pulses reflected off the bottom 
may be used to determine the volume reverberation of bottom material. However, the 
volume reverberation obtained is quite sensitive to inhomogeneities in both small-scale 
bottom topography and sediment composition. 

In fact, even very homogeneous surfaces, such as the air -water interface led to 
slope small fluctuations between successive runs. For this reason, averages of several 
runs were used to determine the volume reverberation level. The use of an average of 
several runs becomes even more important when the sediment grain size increases. 

In the results presented next, this fact is very evident. In Figure 5.1, the decay 
constants obtained for the different sediments show big fluctuations mainly for the 
aggregate and aquarium sand. 

A proportionality between the sediment grain size and the amount of volume 
reverberation was expected as suggested by Clark, Proni, Seem and Tsai (1984), This 
was observed by the averaged final values of the slopes obtained for the four different 


sediments. 


A. FINE SAND 

For the fine sand, eight measurements were made and the individual slopes are 
listed in Table 2. For each slope, the correlation coefficient was also computed. 

Despite the relative stability of the individual slopes, the chain average was 
applied to the 8 measurements and the results are listed in Table 3. 

The running average of all eight data sets is nm = 3.51 x 10° Np/ts and when 
compared to the decay constant of the purely reflected (water-air interface) pulse, @ = 
3.7 7x 10°74 Np/ts, shows a small amount of volume reverberation. This is justified due 
to the small dimensions of the grain size, (0.3 mm) compared with the 8 mm of the 


wavelength, making a a/A = 0.375 which can be considered inside the Rayleigh reginie. 
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TABLE 2 
FINE SAND [INDIVIDUAL SLOPES 


DATA SLOPE x 104 CORRELATION 
SET (Np/is) COEFFICIENT 
DATA41 =6)4 Gk “6, 37S 
DATA42 oe -0. 9162 
DATA43 aes Eig), Sea 
DATA44 -3.58 -0. 9682 
DATA45 -3. 66 -0. 9575 
DATA46 =3EeG ~0. 9526 
DATA47 =3E30 -0. 9485 

| DATA48 23653 -0. 9871 





TABLE 3 
FINE SAND AVERAGE SLOPES 


NO. SETS SLOPE x 107 CORRELATION 
AVERAGED (Np /ts) COEFFICIENT 
2 23570 -0. 9619 | 
3 =2c6 =0.9715 | 
4 aceon oct! 
5 ory, -0.9765 
6 SQWGS 50), Soe | 
7 50 PG: 
g 51 .9757 
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B. AQUARIUM SAND 

A series of ten sets of measurements were run for this sediment (Table 4). The 
individual slopes present a big variation (Figure 5.1). Chain averaging gave the results 
in Table5. 


TABLE 4 
AQUARIUM INDIVIDUAL SLOPES 


DATA SLOPE x 104 CORRELATION 
SET (Np/ps) COEFFICIENT 
DATAS1 = 3800 Tope Ti.ih J 
DATAS2 eo -0. 9424 
DATAS3 -2.68 Bon 7402 
DATAS4 i, aie -0. 6269 

DATASS =2) 43 ops ya’ 
DATAS6 Seely: -0. 9557 

DATA57 =), Of ZOn ST 

| DATAS8 BOG -0. 8495 

| DATAS9 Seley -0. 8564 

DATAG6O =3 Fo), 22377 


The final decay constant Nag = 2.75 x 10° Np/pfis 1s now quite different from 
the reference value of the water-air interface. This difference results from the effect of 
volume reverberation in the sediment and also, as shown in a later measurement, from 
surface scattering due to the relatively big dimensions of its grains compared with the 
wavelength ( a/A = 0.25 ). 
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TABUE 5 | 
AQUARIUM AVERAGED SLOPES 


NO. SETS SLOPE x 10% CORRELATION 
AVERAGED (Np/ps) COEFFICIENT 
2 Eel, ois -0. 9634 
3 =nos -0. 9506 
4 =e -0. 8540 
5 52 OF. ~0. e762 
6 aay 8, S108 
7 -2. 68 -0.9157 
8 -2.68 -0. 8982 
9 eae -0. 9046 
10 =2 75 0), SilZe 





C. AGGREGATE 

Ten measurements were performed using the aggregate with the resulting 
individual slopes listed in Table 6. The variation, in the presented individual slopes 1s 
considerably larger than the variations for the aquarium sand (Figure 5.1), which 1s to 
be expected, since the aggregate shows many inhomogeneities. To make these results 
more understandable the chain average was run again (Table 7). 

The final decay constant Nag = 1.84 x 10° Np/ts, when compared with those 
for the other sediments, shows an even larger amount of reverberation in accordance 


with the proportionality between reverberation level and size of individual scatterers. 


D. GLASS BEADS 
Although the glass beads appear homogeneous in size, shape, and material. the 
individual slopes show considerable variations (Figure 5.1) with the numerical results 


listed in Table 8. The corresponding averaged values are presented in Table 9. 
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TABLE 6 
AGGREGATE INDIVIDUAL SLOPES 





| DATA SLOPE x 104 CORRELATION 
SET (Np/ps) COEFFICIENT 
DATA61 =2eule -0.9764 
DATA62 ago -0. 9385 
DATA63 -2.40 -0. 8392 
DATA64 7) Sail aly 177% 
DATA65 =i), 86 Once 
DATA66 =20S -0. 8901 
DATA67 =2 802 -0. 8187 

| DATA68 =1e55 -0. 9313 
DATA69 snanic Ho) G21 
DATA70 =a, -0. 9612 





The final decay constant for the glass beads is Ny = 3.06 x 10°2 Np/pts which 1s 
consistent with the final slopes of the other sediments, considering their relative grain 


§1ZeS. 


E. OVERALL RESULT 

Relative values of reverberation level for the different materials are shown in 
Figure 5.1 The averaged values of the slopes are marked along the y-axis for the 
different materials. Plotting them against grain size, the corresponding best fitted 
straight line is defined by: 


nN = 3.62 - 0.0035a 


The value of the average grain size of the sediment can then be estimated by : 


a = 10.34 - n/0.0035 


32 





TBA 7 | 
AGGREGATE AVERAGE VALUES 


| 
NO. SETS SLOPE x 10° CORRELATION 
AVERAGED (Np/ts) COEFFICIENT 
2 -2.66 -0. 9476 
3 =2159 -0. 9474 
4 262 -0. 9595 
5 -2.05 -0. 9890 | 
6 2005 -0. 9736 | 
7 Os =Oro730 
8 2700 -0.9756 
9 =), SE -0. 9653 
10 =e “6, SFL 





where s is in Np/fts and a is in mm. This relationship gives only a rough estimate of 
the sediment grain size since, considering the standard deviations in the slopes of the 
extreme data points, the value for a can be estimated within = 2.8 mm, for a slope of 
25x 10 Np/pis. Also this relation was established based on measurements 
performed using a transducer with damping corresponding to a decay constant @ = 
3.77 x 10°2 Np, pis at 180 kHz. Other transducers will have different decay constants 
and further measurements should be performed using different transducers to 


investigate the influence of @ on this result. 
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UUeVoiule & 
GLASS BEADS AVERAGED SLOPES 


NO. SETS SLOPE x 104 CORRELATION 

AVERAGED (Np/us) COEFFICIENT 
2 Bones -0.9794 
3 =e), 1 -0. 9750 
4 =, BE -0. 9646 
5 -2. 86 =09471 
6 a9) ye -0. 9693 
7 -2. 83 -0. 9780 
8 ay ey “05 SORE 
9 =o on -0. 9789 
10 -3.06 -0. 9805 
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Figure 5.1 Spread of individual slopes. 
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Figure 5.2 final averaged slopes vs grain size. 
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VI. CONCLUSIONS 


The purpose of this experiment ts to establish a relationship between acoustic 
volume scattering coefficients and sediment grain sizes. In the experimental design 
only volume reverberation was to be measured, and surface scattering was to be 
considered negligible. Surface scattering was not negligible for most of the materials 
studied, and it was not possible to achieve the initial purpose. 

Nevertheless some conclusions can be drawn: 

1. It was found that sediment grain size is related to the decav constant of the 


tail of the scattered pulse by the equation 
a= 10.34 - n/0.35 (6.1) 


This equation has, however, some limitations since it is applied to the decay constant 
of the specific transducer used and to the frequency of 180 kHz. It can be used, 
however, if the transducer is calibrated to a decay constant @ = -3.77 x 10° Np/ pts. 

2. For wavelengths of the same order of magnitude as the average grain size, 
surface scattering 1s an important component of the reverberation level, even if the 
sediment surface is perfectly levelled and smoothed, and depends also on the sediment’s 
grain size. 

3. The frequency used and materials selected for the experiment correspond to a 
wide range of ka values, from 0.12 for the fine sand to 4.16 for the aggregate. The 
decaying slope monotonically decreases with increasing grain size and does not show 
any strange behavior in the region of ka = lI. 

4. Inhomogeneities in sediment composition are not responsible for the variation 
in the decaying slopes obtained, at least not totally, since measurements made on 
homogeneous material like glass beads show the same spread in results as the aquarium 
sand of similar size. 

Further study of the volume reverberation from ocean sediments is 
recommended, as is the use of different combinations of frequency and grain sizes, 
keeping in mind that wavelengths of the same order of magnitude as the grain size 


should not be used so that surface scattering may be neglected. 
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A theoretical search must be done to Investigate the influence of transducer 


characteristics on the decay constant, to verify the extent of the validity of Equation 
6018 


APPENDIX 
DATA COLLECTING AND PROCESSING SEQUENCE 


In this appendix a complete sequence of the data collecting and processing for 
the fine sand is presented. Eight data sets (DATA#) of 50 pulses each were collected 
for this sediment and averaged using the program THESIS3. These data sets (files 
DATA4I to 48) are graphically represented in Figures 1 to 8 where the squared 


voltages are plotted against time (program PLOT). 
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Figure 1 Pulse waveform of DATA4I. 


The decaving slopes for these pulses are then determined using the program 
MAXIMUM that selects the points of local maximum voltages. A list of these points 
on the pulse tail for the first data set (DATA4IL)is presented next. 
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LISTING OF LOCAL MAXIMA FOR DATA4I 
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The corresponding waveform envelope is plotted in Figure 9. 


connecting the points of maximum voltage. 
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It 1s obtained by 
ies best titted) straight line and 


corresponding slope are then determined (program PLOT3). 


The individual slopes show big spread (sce Chapter five), and they were applied 


in the averaging process. The eight data sets were averaged as follows: 


DATA48 
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Figure 2. Pulse waveform of DATA42. 
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Figure 3 Pulse waveform of DATA4S3. 
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Figure 4 Pulse waveform of DATA44. 
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Figure 5 Pulse waveform of DATA4S. 
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Figure 6 Pulse waveform of DATA4G6. 
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Figure 7 Pulse waveform of DATA47. 
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Figure § Pulse waveform of DATA4S8. 


DAT42 1s the average of DATA4I and DATA42. From DAT42 to DAT48 the 
Variations presented in the individual waveforms are filtered because the number of 
terms involved in the average increases. DAT48, the result of the averaging, may be 
considered the representative waveform for the fine sand. 

Next, programs MAXIMUM and PLOT3 are apphied to the averaged DAT# sets 
which corresponding plots are presented in Figures 10 to 17. 

The envelope from DATA4I to DAT48 gets smoother and the value of the decaying 
slope stabilizes along the averaging process although the differences are small since the 


individual waveforms are not very different for this sedirnent. 
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Figure 10 DAT42. 
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Figure 11 DAT43. 
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Figure 12 DAT44. 
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Figure 13. DAT45. 
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Figure 14 DAT46. 
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Figure 15) DAT4?7. 
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